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Abstract 

We investigate the dynamics of a particle executing a general Continuous Time Random Walk 



(CTRW) in three dimensions under the influence of arbitrary time- varying external fields. Contrary 
g to the general approach in recent works, our method invokes neither the Fractional Fokker-Planck 

■d^ equation (FFPE) nor the Stochastic Langevin Equation (SLE). Rather, we use rigorous probability 

c/3 arguments to derive the general expression for moments of all orders of the position probability 

I ^..^.^^^„^.„,,^.„^.,^,„^^^, 

1 Closed form expression for the position probability density is derived for the memoryless condition. 

^ For the special case of CTRW on a one-dimensional lattice with nearest neighbour jumps, our 

I ^1 equations confirm the phenomena of "death of linear response" and "field-induced dispersion" for 

^ sub-diflFusion pointed out in [I. M. Sokolov and J. Klafter, Phys. Rev. Lett. 97, 140602 (2006)]. 

However, our analysis produces additional terms in the expressions for higher moments, which 
have non-trivial consequences. We show that the disappearance of these terms result from the 
approximation involved in taking the continuum limit to derive the generalized Fokker-Planck 
equation. This establishes the incompleteness of the FFPE formulation, especially in predicting 
J> the higher moments. We also discuss how different predictions of the model alter if we allow jumps 

^ beyond nearest neighbours and possible circumstances where this becomes relevant. 
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I. INTRODUCTION 



In the past few decades, an increasing number of natural phenomena in diverse branches 
of physics, biology and the related fields of science have been shown to exhibit anomalous 
diffusion. Examples of processes displaying subdiffusive dynamics include gel electrophoresis 
of polymers (e.g., DNA), transport in amorphous semiconductors [HE], nuclear magnetic 
resonance, motion of proteins through cell membranes, diffusion in percolative and porous 
systems, transport in fractal geometries, contaminants transport by underground water and 
many others (see P] for interesting reviews). The Continuous Time Random Walk (CTRW) 
model introduced by Montroll and Weiss [U-E], with power-law waiting time distributions 
lacking first moment, provides a powerful mathematical tool to describe sub-diffusion [U 
[7j. The non-stationary nature of these (non-Markovian) subdiffusive CTRWs leads to the 
effects of aging in various systems pHTO]. In the absence of time-dependent external fields, 
the CTRW process can be subordinated to the simple random walk, which enables one to 
describe it within the framework of Fractional Fokker-Planck Equations (FFPEs) |31 ttU- 
The FFPE can also be derived from the analysis of a (ordinary) Langevin equation in 
subordinated random operational time [I5l [16] . 

However, the generalization of the above scheme to situations where an external time- 
varying field is present is non-trivial because the force physically changes in real time which 
cannot be transformed to the random operational time, thus creating difficulties in subordi- 
nation [T7l[IS]- It has been shown that the the FFPE for the time-independent case, when 
generalized ad hoc to a time-dependent force [HI ^U\, does not correspond to physical real- 
ity [18]. Modified FFPEs for time- dependent forces were derived in Refs. [iHl [19] but later 
more rigorous treatments revealed that they were valid only for some particular deterministic 
time-dependences [201 [21]. Recently, a subordinated stochastic Langevin equation [22j and 
an FFPE [23] have been proposed for modelling sub-diffusion in space-and-time-dependent 
force fields. 

In spite of these advances, the theoretical investigation of sub-diffusion in a time- 
dependent force field still presents a challenge that is far from trivial. Some of the diffi- 
culties and ambiguities were pointed out in Ref. [21]. First, it is known that there is no 
unique physical mechanism behind the occurrence of sub-diffusion in condensed media [23] . 
Second, in the CTRW model, successive waiting times are assumed to be mutually uncor- 
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related, which may not be true in practice [25]. Third, as the time- varying field will alter 
the rates for moving in different directions when escaping from a trap location, the waiting 
time distribution will generally be affected by the field [2B]. However, due to our lack of 
knowledge of a concrete mechanism, it is not clear to what extent this infiuence will change 
the predictions based on a decoupled model where the waiting time distribution is assumed 
to be independent of the external perturbation [27]. 

The limitation we want to stress in this paper is that even for an uncorrelated decoupled 
CTRW model, the FFPE formulation is incomplete [21], as no (single event) non-Markovian 
master equation can fully characterize the underlying (multitime event) non-Markovian 
stochastic process which lacks the factorization property [28j. Therefore, neither a non- 
Markovian FFPE nor its solution can display all the features of the non-Markovian CTRW 
[29] . Further, the general strategy to form an FFPE adopted in all relevant works (see, for 
example [19]) is to first write down the generalized master equation (GME) for a CTRW 
on a one-dimensional lattice with jumps restricted to nearest neighbours, and then take 
its continuum limit. This method is not rigorous and as we show in this article, the pre- 
dictions of the FFPE regarding the higher moments of the position probability density of 
the random walker increasingly differ from those of the GME. Also, if the CTRW is not 
confined in a single dimension, there are several non-unique ways of taking the continuum 
limit corresponding to different lattices. 

In this paper, we treat the problem of purely time-dependent external fields. Contrary to 
the customary approach, we start with a general CTRW in a three-dimensional continuum 
and derive expressions for the moments of all orders of the position probability density 
p{r,t) of the walker for arbitrary waiting time density iplt) and jump probability density 
0(r|t). Later on, we choose a particular form of 0(r|t) to recover the special case of the 
one-dimensional lattice. We find that in this case, our equations do predict the same mean 
position and standard deviation as those in Refs. [TSH2T] . thus confirming the phenomena of 
"field-induced dispersion" and "death of linear response" for systems showing sub-diffusion. 
However, for the higher moments, our analysis produces extra terms in addition to those 
in the aforesaid references, which have non-trivial consequences. We point out that the 
disappearance of these additional terms is an artifact of the approximation involved in 
taking the continuum limit to derive the generalized FFPE. These terms are recovered if 
one calculates moments directly from the GME. Our method is based on direct probability 
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arguments [30] and some well-known properties of Integral transforms and does not invoke 
the fractional Fokker-Planck equation. This is a key departure from other recent works on 
the subject. It presents a rigorous treatment of the general CTRW in three-dimensions. 
Consequently, the results are valid for any time-dependence of the external field within 
the assumptions of uncorrelated decoupled CTRW. We find the interesting result that the 
nth moment of p{r,t) depends on the mth moment of (f){r\t) if and only if m < n. For 
sub-diffusive waiting time densities, tplt) ~ t~^^" with < a < 1, we show how the normal 
approximation breaks down more and more as one lowers a towards zero. We describe a one- 
parameter family of probability density functions which can well-approximate the position 
probability density of the random walker for large values of t. For the special case where 
the random walker is memoryless, i.e., where the probability of making a jump does not 
depend on the time at which the previous jump was made, the process becomes Markovian. 
In this instance, we are able to derive p{f, t) in closed form and show how its moments are 
simply related to the moments of (f){r\t). We conclude with a discussion on the significance 
of choices of for which the jumps are not restricted to nearest neighbours and a brief 

note on some open problems in the subject. 

II. FORMULATION OF THE PROBLEM 

Let us first describe the model of the three-dimensional random walk adopted in this 
paper. The random walker, on arriving at a position p at time r stays there for a sojourn 
time t\ after which it jumps to a position r at time t = t + t' . The waiting time t' 
is distributed with the probability density function ^p{t'). The jumps are assumed to be 
instantaneous on the time scale of typical waiting times and changes of external parameters. 
A physical justification of this assumption may be found in the case of sub- diffusion, where 
long waiting times are often dominant. The probability that a jump occurring at time t 
results in a displacement between r' and r' -|- dr' is given by (f){r'\t)d^r' . We seek to find the 
probability density p{f, t) that the random walker is found at position r at time t. Also, we 
choose our references in such a manner that at t = 0, the walker has just arrived at r = 0. 

Let us define as the probability that the time interval between two successive jumps 
is greater than t. Clearly, 




(1) 
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We shall require another intermediate probability density f{r,t) defined by, 

f{r,t)d^fdt = probability that the walker arrives between f^and f+df 

within the time interval t to t + dt after one or more jumps. 

Then the desired probability density p{r, t) can be written as a sum of two parts, 

p{f,t) = 6{f)^{t) + q{f,t) , (2) 

where, q{r,t) = / f{f,T)'^{t — T)dT. (3) 
Jo 

The first term involving the delta function expresses the possibility that no jump has 
taken place since t = 0, in which case, the random walker is bound to remain at = 0. 
Conversely, the probability density of being at r at time t after executing at least one jump 
since t = is given by q{f, t), which can be written as a convolution integral as in (g. Here 
r denotes the time of occurrence of the last jump which took the walker to f^and the second 
factor "^(t — t) ensures that no jump has taken place since then. From ^ and ([s]), it is 
evident that our problem reduces to the determination of the intermediate density f{f,t). 

This intermediate density can again be broken into two parts. It is possible that the 
only jump took place at time t, taking the walker to r, for which the probability density is 
given by ip{t)(j){f\t). Alternatively, the walker reached some position p at some time r lying 
between and t, and was transferred to r in the last jump at time t. This leads us to the 
following integral equation for f{f,t). 

/oo 
/ fip,T)^Pit-T)(l)if-p\t)drd'p. (4) 
oo Jo 

This equation can be simplified by taking the Fourier transform with respect to r and 
using the convolution theorem, 

m t) = mt)m + mt) f m r)Ht - r)dT . (5) 

For the special case where 4>{r\t) is time-independent, it is possible to solve this equation 
simply by taking the Laplace transform with respect to t, which then leads directly to the 
Montroll- Weiss equation [4J. However, for a non-trivial time-dependence in (j){f\t), no such 
straightforward method seems to exist. In the following, we first assume a specific form of 
the waiting time density. After discussing its implications, we proceed to the general case. 
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III. RESULTS FOR MEMORYLESS RANDOM WALK 



-t/T 



Using this result in ([T]), we obtain the waiting time density 



In a memoryless random walk, the probability of a jump in any infinitesimal time interval 
dt is a constant and given by dt/T, where T is the average time interval between two 
successive jumps. Then the probability that no jump occurs in an interval of duration t is 
given by, 

(6) 



(7) 

(8) 
(9) 

(10) 
(11) 



dt ' ' T 
Substituting the above in ^ and ([s]), we get. 



qik,t) 

mt) 



Jo 

1+ [ f{k,T)e^/^dT 
Jo 



From these two equations, one can readily derive the following. 



qik,t) 



-t/T 



T 



(p{k\T) [l + q{k, T)e^/^]dT . 



Differentiation with respect to t yields. 



^^q{k,t) + ^{l-mt))q{k,t) 



T 



e-*/^0(fc|t) . 



This can easily be solved using standard methods. After substituting the boundary 
condition q{k,0) = 0, which follows from (pj), we obtain. 



q{k,t) 



-t/T 



exp I ^ / 0(A;|r)rfr) - 1 



(12) 



To obtain p{k, t), we take the Fourier transform of (|2j) and substitute the above expression 
for q{k,t). This gives. 



p{k, t) = e ^^'^ exp 



T 



(j){k\T)di 



(13) 



This formula can be verified in some simple cases. First, for very small values of t, this 
can be expanded as, 

i^(^,t)=(l-^)(l + ^0(^|O))+O(t^) 
1-^).1 + ^0(^|O) + O(t^). 



Taking the inverse Fourier transform, one gets, 

p(f, t) -- 



(14) 



no jump one jump 

A little reflection reveals that the first (or second) term on the right hand side is nothing 
but the contribution of the possibility that no jump (or a single jump) occurs in time t. 

Consider another simple case where every jump results in a constant displacement a of 
the walker, i.e., (f){r\t) = 6{f — a). Then 0(/i;|t) = e"**^'^ and (13) gives, after expansion of 
the exponential. 



n=0 

oo 



Therefore, p{r,t) = / ^ — mr, ^{^ ~ ^^) 



n=0 



n! T" 



(15) 



Now, the probability that exactly n jumps take place in duration t is given by, 
T^nit) = / ... / / V^(ri)V^(r2-ri)...V^(r„-r„_i)^(t-r„)rfrirfr2...rfr„_i(ir„ 



JO 



^0 



(16) 



t/T ft pTn 



^0 



T3 rT2 



^0 



dTidT2 . . . dTn-ldTn 



n! T" 



Comparing this with (15), we see that p{r,t) can be written as, 

oo 

p{r,t) = ^Vn{t)5{f - no) , 



n=0 



exactly what it should be. 



Thus, for the memoryless random walk, given the jump probability density 0(r|t), (13) 



can be employed to calculate p{k,t). Then, in principle, p{r,t) can be obtained via the 
inverse Fourier transform. However, because (f){k\t) is present as an exponent, even for 
simple choices of (j){f\t), the inverse Fourier transform may not be available in closed form. 
Nevertheless, moments of all orders of the distribution p{f, t) can be calculated from p{k, t) 
itself and these moments are what we can measure experimentally. 
From the expression of the Fourier transform 



p{k,t) 



one can write the nth moment of p{f, t) as, 



< x" > 



(9" _ - 
x^p{f,t)dx = i^—p{k,t) 



k=0 



(17) 



Inserting the expression of p{k,t) for memoryless random walk (see (13)), one finds 



< x" > 



k=0 



where, g{k,t) 



T 



(j){k\T)dT 



(18) 
(19) 



The nth derivative in (18) can be written as a sum of terms involving mth derivatives of 



g{k,t) with m <n (assuming both m and n to be nonnegative integers), as in the following, 

n = 1 : e^g' 

n = 2 -.e^g'^ + g") 

n = 3 ■.e^{g" + ^g'g" + g"') 

n = 4 : e^{g'' + Qg'^g" + Ag' g'" + ?,g"^ + g"") 

: etc. . 



(20) 



In these expressions, g' denotes the derivative of g{k,t) with respect to k^. Adopting 



similar conventions, we can write, from (19) 



7HI 



IA;=0 



1 1 



< >s dr where, < x" >, 







x"0(f|t)(ia; . 



(21) 



In deriving the second equality, we have made use of the relation expressed by (17). Also, 

ft 



T 
1 
f 



0(O|r)c/r 







t POO 



{r\t)d^fdT 



t 

f 



(22) 



'0 J-00 

The last equality follows from the fact that 0(r|t) denotes a probability density and 



hence (f){r\t)d^r must be equal to 1. Equations (18) to (22) enable us to write down the 



following compact relations which illustrates very clearly how the moments of p(r*, t) depend 
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on those of the jump probabihty density. 



< X > = < X >^ 

< x^ > = < X >^ + < a;^ >0 Thus, standard deviation Ax = \/ < x"^ ><, 

< > = <1F>0^ + 3<ir><^< x^ >0 + < x^ >0 (23) 

< x^ > = < X + 6< X >(/,^< x^ >0 + 4< X >0< x^ >0 + 3< x^ ><^ + < x"^ ><^ 

: and so on. Here, a = ^ /J a{T)dT. 

Similar equations exist for y and z as welL These lead us to the following expressions for 
the mean position and standard deviation in three dimensions. 

<f> = <x>i+<y>j+<z>k = <f>^. (24) 



Ar = \/<r^>—\<f> 



< x^ > + < y^ > + < z^ > -\< r >^ 



< r2 >^ . (25) 



The relations in ( 23 ) can be more meaningfully expressed in terms of the central moments 



fin defined as, fin =< {x— < x >)" >, as follows, 

/ii = < X >0 , fi2 = < x'^ >^ , /i3 = < x^ >^ , /i4 = < x"^ + 3< x2 >0 etc. (26) 

The significance of the second term in the expression for fx^ can be understood in terms 
of the kurtosis k = fii/ fi\ — 3. The kurtosis is a measure of the peakedness of a distribution. 
Higher k values result in a sharper peak with heavier tails. The normal distribution possesses 



K = 0. From (26), one finds that, n = < x^ ><t>/< ^"^ ></> ? which is always positive. Thus, 



the position probability for the memoryless random walk will always be more sharply peaked 
than the normal distribution. However, if < x'^ >^ and < x^ >^ are time-independent, then, 
K will approach zero for large values of t. 

Also, if the jump probability is symmetric such that < x" >0 vanishes for all odd values 



of n, then, from (18) and (20), it is evident that all odd moments of p{r,t) will vanish as 



well. Therefore, the position probability density will be symmetric, as expected. 
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IV. RESULTS FOR GENERAL RANDOM WALK 



Here we proceed to the more challenging task of deriving the moments for a general wait- 
ing time density. We start by expressing the moments of p{f, t) in terms of the intermediate 
density f{f,t) and waiting time density ^p{t). Using ^ and (|3]), one can write, 



< a; > (t) = / xp{f^t)d^f 



— oo 
oo 

xq{r, t)d^r 



iq'{0,t) using (17) . 



= if'{0,t) * ^(t) where * denotes convolution in time domain. (27) 

Taking the Laplace transform of both sides, one obtains, < x >{s) = if'{0, s)^(s). But, 
according to the definition of ^(t) in (ll|, ^(s) = (1 — ijj(s))/s. Hence, we get. 



<x>{s) = zf{0,s)(^-^^] . (28) 



Similarly, one can show that, 

< >{s) = -f"iO, s) (^^^) , (29) 
and, in general, < af" >(s) = ^/("^(O, s) [ i^^l^ ) , n = 1, 2, 3, ... . (30) 



Therefore, in order to determine the moments of p(f,t), we must first find out /'•"'' (0, s) 
from the integral equation ([s]). For the purpose of simplifying mathematical notations, we 
define a linear operator A by stating its action on an arbitrary function g{k,t) as, 

A[g{k,t)] = mt){9{k,t)*ij{t)} . (31) 

Equation (|5| can then be written in a compact form, 

f{k,t)-A[f{k,t)] = mt)m- (32) 

Here we note a property of the operator A. Consider any two functions h{k, t) and g{k, t) 



related by h{k,t) = A[g{k,t)]. Using (31), one can write 



h\k,t) = '^^hik,t) + A[g'ik,t)]. (33) 

nk\t) 
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Taking derivatives of both sides in (32) and applying the above property, we obtain, 



(j){k\t) 

or, f{k, t) - Aifik, t)] = t) . 

(f){k\t) 



(34) 



Using (17) and the fact that (j){r\t) is a probabihty density, one can write, 

POO 

(j)(0\t) = / (j){f\t)d^f= 1 , 



and, (p'iOlt) = - < X >s (t) . 
i 



(35) 
(36) 



Again, it follows from (31) that 



A[g{k,t)]\^^^ = g{0,t)*m ■ 



(37) 



Application of these results on (|32|) and (|34|) yields, 

M^(s) (say) 



and, f(0,s) 



Substitution of ( 39 ) into ( 28 ) gives 



< X >{s) 



C[<x>^ {t)M4t)] 



(38) 
(39) 



(40) 



Taking inverse Laplace transform, we get the general expression of the mean position. 



< X > (t) = < x>^{t) M,^{t) dr . 



(41) 



Combining analogous expressions for < y > and < z > and making use of (38), we can 
thus write, 

' 1 ir . (42) 



<f>{t)= < f>^ (r) C-^ 
Jo 



This can be explained as follows. From (16), we find that the probability of exactly n 



jumps taking place in duration t can be written as. 



^ V ' 

(n + 1) convolutions 



(43) 
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Therefore, the mean number of jumps occurring in time t is given by, 



< N{t) > = C-^ 



n=0 -I 

f:n«.(.))"(i^ 



n=0 



(44) 



Thus, on average, the number of steps taken by the random walker in an infinitesimal 
time interval dt is, 



d < N{t) > 
dt 



dt = c-^ 



dt . 



(45) 



Now we see that the integrand in (42) is the product of the average number of jumps in 



duration dr and the average displacement of the walker in each such jump. So, it represents 
the average displacement in the interval dr, which when integrated, gives us the mean 
position of the random walker after time t, as expected. In fact, this reasoning was used to 
derive the linear response of systems to time-dependent external fields [9]. 

For the special case of the memoryless random walk, ilj{t) = (1/T) exp(— t/T). Hence, 



ij{s) = l/{Ts + 1), giving M^(s) = 1/Ts and M^{t) = 1/T. Putting this into (|4l|), we get, 
< X > (t) = (1/T) < X >^ (r) dr, in agreement with what we had derived in (23). 



Next we proceed to calculate the higher moments. Keeping analogy with (32) and (34) 
let us write, 



f-\k,t)-A[p'^\k,t)] = Xn{k,t) 



f{") 



n = 1, 2, 3, - 



Then (30) and (37) together imply. 



< X" >{s) 



-Xn(0, S) 



or, -<x^>{t) = z"xn(0,t) 



(46) 

(47) 
(48) 



Xn{k,t) can be obtained by taking the n-th derivative of (32) and substituting the result 



in (46), 



Qn 



Xnik,t) = — (0(%){Mt)*V^(t)}) -0(fc|t){/(")(^,t)*^(t)} + 0('^)(fc|t)V^(t) . (49) 



12 



Expanding the first term by Leibnitz's formula, we get, 



9=1 



n-1 



q=l 

Tlie term outside summation can be simplified for /c = by noting that. 



(50) 
(51) 



C[iP{t) + f{0,t) * tp{t)] = ^(s)(l + M4s)) = M^s) [using (38)] 



(52) 



This result, along with (17), (30) and (48), yield the following equation giving the n-th 



moment in terms of the lower moments, (n = 2, 3, 4, . . . ) 

n-l 



d 

di 



< > (t) = < >^ {t)M4t) + < >0 W-^"' b < >is)M4s)] (53) 

q=l ^^-^ 

n~l / s , 

= <x^>^ it)M4t) + 5^ r j < x'^ >0 {t){M4t) * - < x^-" > (t)} . 

(54) 

For the memory less random walk, we noted that < x" > depended on < >^ if and 
only if m < n. According to the equation above, this fact remains true in the general case. 



It can be readily verified that (54) reproduces the results for the memoryless random walk 



for the choice iplt) = (1/T) exp(— t/T). As we found earlier, for this waiting time density, 



M^{t) = 1/T. Substituting this in (54), we obtain, 

> (t) = < >^ (t) < x-^ > (t) 

q=l 



< X 



(55) 



with the understanding that < x° > (t) = 1. That this equation is indeed equivalent to (18) 



can be confirmed as follows. Differentiating (18) with respect to t and using (19), one can 
write, 



d 
dt 



1 1 (9" r 

< > (t) = - - < > (t) + -z"e~*/^^ 



^{k,t) 



k=0 



1 

f 



q=0 



dkx 



k=0 



q=l 



(56) 



13 



which is identical to (55). 



Equation (54) and (41) constitute the complete solution for the moments < > (n G 



Z+) of p{r, t). These equations, along with their counterparts for y and z, completely specify 
the random walker's position probability density p{r,t). 



As a special case, let us investigate the standard deviation Aa;. For n = 2, (54) becomes, 



after substituting from (41), 

^ < > (t) = < a;' >^ (t)M^(t) + 2<x>^{t) {M^t) *{<x>^ (t)M^(t))} . (57) 
Hence, 

^ < {Axf >{t)=j^<x^>{t)-2<x>{t)j^<x> it) 



< x^ >^ {t)M4t) + 2<x>^{t) {M^t) *{<x>^ (t)M^(t))} 



-M^t) / <x>^{t)M4t) dr 



(58) 



Thus, the standard deviation is given by. 



< x^ >^ {r)M^{T) dT+ 2<x>^ (ri) {M^n) * {< x >^ (ri)M^(ri))} 



Ax 



- M^(ri) / <x>^ (r2)M^(r2) dT2 dn . (59) 
Jo 

If the random walk is memoryless, M^{t) = 1/T, a constant. As a result, the second part 
in the above equation vanishes and we get. 



(^^)memoryless \J T Jq ^'^ ^<t> '^^ 



(60) 



conforming to what we found earlier (see (23)). 



V. DISCUSSION OF RESULTS 

A. A Case Study : Failure of generalized Fokker-Planck equation in predicting 
higher moments 

A few years ago, Sokolov and Klafter[T9] analyzed a similar problem where the random 
walk is one- dimensional and the walker takes steps of unit length in either direction. This 
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is seen to be a special case of our problem where the jump probability is given by, 

<f>im = Q + ^/(t))5(r - 5) + - f + a) • (61) 

where a is a unit vector along positive x direction. Consequently, the moments are given by, 

< ><^ = /i/(t) if n is odd , (62) 
= 1 if n is even . (63) 



When these are substituted, Eqs. (41) and (57) indeed reproduce Eqs. (12) and (13) in 
Ref. [19]. Thus, the expressions we deduced for the mean position and standard deviation of 
the random walker agree with those in Ref. [19] found via a different method. Consequently, 
the comments made there regarding the phenomena of "field-induced dispersion" and "death 
of linear response" in systems showing sub-diffusion apply here also. However, the higher 



moments obey the following equation, obtained after substituting (62) and (63) into (53) 



|<a:">(t)= ('')^<x''-'>it)+fif{t) (""Y < ^""-^ > it) . (64) 

9=2,4,6,... (j=l,3,5,... 

where $ < a;"^'^ > (t) = C^^[s < a;"^'? >(s)M^(s)] coincides with the integro-differential 
operator used in Ref. [19j. The general equation for moments derived in Ref. [19] (see 
equation (11)) from a generalized Fokker-Planck equation (FPE) contains only the first two 



terms corresponding to g = 1 and q = 2 in (64). Similar equations were derived in Refs. 



[TS| f2U\ f2T\ . The two equations match only if we substitute < x"^ >^ = 0Vra>3 (which 
is not possible for any choice of (j){f\t)). However, we note that the loss of the additional 
terms is an artifact of the approximation involved in taking the continuum limit of the 
Generalized Master Equation (GME) to derive the generalized FPE. As illustrated in Ref. 
[19j . application of probability conservation arguments yields the following GME, 



Pkit) = Q + ^/(t))$P.-i(t) + - ^/(t))$P.+i(t) - hkit) , (65) 
where Pk{t) denotes the probability that the random walker is at site k at time t. The n-th 



moment can be calculated as, < x" > (t) = J2'k=-oo^"'Pki't)- Thus, multiplying (65) by k' 
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and summing over k, we get, 



dt 



^ < > (t) = - $ < x'^ > (t) + Q + ^f2{k + l)>,(t)] 

^ ^ k=-oo 

1 \ °° 

' fc=— oo 

- $ < X" > (t) + Q + ^mYijz Q < ^""^ > w] 



5^ Q^<-"-^>w+/^/w E il)^ 

=2,4,6,... (7=1,3,5,... 



< > (t) , (66) 



which is identical to (64). This clearly shows that the approximations involved in taking the 
continuum limit lead to the disappearance of these additional terms. Although this does 
not affect the mean and standard deviation, the higher moments are increasingly affected. 

The significance of the additional terms in our analysis may be investigated by examining 
how they modify the asymmetry and peakedness of the position probability density p{f,t). 
For example, the 3rd central moment, given by /is = < (x— < x >)^ >, is a measure of the 
lopsidedness of a distribution. This can be expanded to give, fi^ = < > —3 < >< x > 



+2 < X > . Thus, the extra term in (64) for n = 3 increases the 3rd central moment by an 
/o 



amount A/is = 1^ Jq f{T)M^{T)dT. A quick reflection reveals that this is just the expression 



for the mean position < x >. The skewness 7 is defined by, 7 = (^^/{Ax)^. Consequently, 
the additional term results in an increment of 7 given by, 

A7 = — — . (67) 

A positive skewness implies that the probability density function has a longer tail to the 
right of the mean so that the bulk of the values lie to its left and vice-versa. Similarly, 
a measure of the peakedness of a probability distribution is the kurtosis defined as, k = 
Hiji^AxY — 3, where /i4 is the 4th central moment, given by, /i4 =< (x— < x >Y >. Higher 
value of K indicates that the distribution has a more acute peak around the mean and 
fatter tails. Therefore, more of the variance is the result of infrequent extreme deviations. 



as opposed to frequent but modest deviations. The additional terms in (64) increases the 
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kurtosis by, 



2 < a;2 > -4 < a; >2 - /q M^{T)dT 



(68) 



(Ax)4 

Let us examine the order of magnitudes of these changes for the particular case of os- 
cillating external field f{t) = /o sin cut and subdiffusive waiting-time density (ipit) ~ t~^^" 
with < a < 1) investigated in Ref. [TDj. By using the shift and convolution properties 
of Laplace transforms and taking limit as s — )■ 0, one can show that for large values of t, 



< a; > ^ /i/o lm{M^{—iuj)} , a constant and Afi^ increases as (see (A3) and (B9)). Again, 



From (A3) and (A7), we see that for large values of t, the variance behaves as. 



(Ax)2~ (l+/iVoM^l}) / M4T)dT -^VoM-2l^2}-/iVo (M^l})' 

Jo 



(69) 



which increases as t". Substitution of these results in (67) and (68) yields A7 ~ t 



and A/t ~ t However, since we are concerned in how these additional terms modify the 
results in Ref. [19j, the quantities of interest are not the absolute increments, rather the 



relative increments, A7/7 and Ak/k. As we show in the Appendices (see (A12) and (Bll)), 
for large values of t, 

A7 Imjzil 

7 " 6Re{^i}Im{^i} - lfi^f^lm{z^z,{z, - z^)} + 3fi^f^lm{z,}Re{z,Z2} + 2fi^fi{lm{z^}Y 

(70) 

Ak (1 + V/oRe{^i}) /o M^{T)dT + C 



and, 



^ 6(1+ /i2/oRe{zi})' /o M^t) * M4T)dT + A £ Mv,(r)dr + B 



(71) 



where zi = M^{—ioj), Z2 = M^,{—i2ijj), z^ = M^p^—iSu), Z4 = M^(— ■j4(X') and A, B, C are 
constants related to Zi {i = 1,2,3,4). For subdiffusive waiting time density iplt) ~ 
with < a < 1, f^M^{T)dT ~ t° and /J M^(r) * M,^{T)dT ~ t^a^ Therefore, for large 
values of t, A'-f/'-f approaches a constant and Ak/k ~ t^". Thus, in terms of their effect on 
peakedness of the position probability density, the additional terms become more significant 
for smaller a. 

It is also of interest to find out what value 7 and n assume for large values of t. This can 
be useful to determine to what extent the position probability of the random walker differs 
from the normal distribution, which has both 7 = and k = 0. First, we have already noted 
that A7 ~ f-3a/2 g^j^^ A7/7 approaches a constant. Therefore, we conclude that 7 must 
decay towards zero as However, both Ak/k and Ak decrease as Thus, for large 
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t, the kurtosis approaches a constant. From (BIO) and (69), one finds that this constant is 
given by, 

(72) 



K ~ — z ?, 6 



{J^M4T)dTY 

Let us take M^(t) = corresponding to subdiffusive waiting time densities. Then, 

M^{s) = r(a)/s". Therefore, M^t) * M^{t) = C-^[{M^{s)f] = {{T{a)y /T{2a))t'^''-\ 



Substituting these into (72), we obtain, 

K ~ 3 



(r(«))^ ; 



(73) 




a 



FIG. 1: Variation of kurtosis for large t with a. 

Figure [T] plots this hmiting value of k as a function of the exponent a. For a = 1, 
K,{t — i- oo) = 0. This is expected because we know that a = 1 reproduces normal diffusion 
for which the position probability of the particle approaches a Gaussian distribution for large 
t. The figure depicts that as we lower a towards zero, K,{t oo) rises monotonically from 
towards 3. This means that p{f,t — oo) becomes more and more sharply peaked with 
fatter tails. The probability of infrequent large deviations from the mean increases while 
that of frequent modestly sized deviations decreases. This clearly shows how the normal 
approximation breaks down more and more as one enters deeper into the subdiffusive regime. 
The figure also provides us a good picture of the position probability density in the large 
time limit in terms of well-known distributions for certain values of a. For example, the 
logistic distribution has zero skewness and a kurtosis of 1.2. From Figure [T| we see that 
K{t —7- oo) = 1.2 for a ~ 0.64. Thus, we can say that for a ~ 0.64, the pdf of the logistic 



distribution with mean and variance given by (A3) and (69) respectively provides a good 
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approximation to p{f,t — )■ oo). It describes accurately the first four central moments and 
since the higher moments generally measure properties further and further out in the tails, 
we can get a fair idea of its overall shape, especially near the mean. Similarly, ]3(r, t — ?■ oo) 
is well approximated by the pdf of the Laplace distribution for a — )• 0. 
Now consider the probability density function 



/3 



26i//3r(i//3) 

It has the following characterizations. 



exp 



\x — jJ,\ 



(74) 



< X > 



7 



/i 
. 



(Ax) 



r(3//3) 



K 



r(i//3) 
r(5//3)r(i//3) 



3 . 



(75) 
(76) 



(r(3//3))2 

K continuously increases from to 3 as the parameter /3 decreases from 2 to 1. Therefore, 
for any value of a G (0, 1), there exists a unique value of /3 G (1, 2), such that fp{x\fi, h) has 
the same kurtosis as p(r, t — )■ oo). This is illustrated in Figure |2} 




a 



FIG. 2: /3 vs a: The curve is well approximated by /3 = a + 1. 



Thus, for suitable choices of ii and 6 in (74), we can get a rough picture oi p{f,t — )■ oo) 
for any a G (0, 1). 



B. Other Issues: Choice of the jump probability density 



The particular form of the jump probability density given by (61), assumed in recent 



works, implies that < x" is independent of the external field for all even values of n. 
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This may not be true in all circumstances, e.g., when the external field is temperature. An 
increase in temperature is not expected to make a jump towards right more probable than 
a jump towards left. Rather, the increased thermal energy of the particle should allow it 
to hop longer distances, thus increasing < >^. This can arise, for example, if a rise in 
temperature increases the width of (pirlt) which has the shape of a Gaussian distribution 



centered at the origin. In such a circumstance, < x ></>= 0. Therefore, from (41), < x >= 0. 
However, < x"^ >^ {t) = g{T{t)), where g is an increasing function of temperature T. Hence, 
from (57), < x^ > (t) = g{T{T))M^{T)dT. If, for instance, g is a. linear function of 
T which in turn rises linearly with time, then the standard deviation of the particle will 
increase as for subdiffusive waiting time densities, although its average displacement 

will be zero. 

Further, the field-induced component of dispersion in the case of periodic zero-mean 
external field will, in general, depend on 0(r|t). For example, if 4>{'r\t) = (l/v^27r) exp(— (x — 



/i/(t))^/2), then, < x >^= ^f{t) and < x^ >^ (t) = 1 + fi^{f{t)y. Using (57), one finds that 
the dependence of < (t) on f(t) enhances the field-induced dispersion by an amount 

^'^{f{T)YM.^{T)dT. For sinusoidally oscillating field /(t), this contribution grows as 
for large values of t, which is at the same scale as the original estimation in Ref. [19] based 
on a nearest-neighbour jump model. Thus, relaxing the assumption of nearest neighbour 
jumps can significantly alter the predictions of the model. 



VI. CONCLUDING REMARKS 



In this paper, we derived rigorous results for a general CTRW in three dimensions under 
time-dependent external forces. The method is based solely on probability arguments and 
does not invoke either the FFPE or the SLE. In this sense, the approach is more direct 
and illuminating. In the special case of CTRW on a one-dimensional lattice with nearest 
neighbour jumps, our equations confirm the phenomena of "death of linear response" and 
"field-induced dispersion" for zero-mean periodic forcing with subdiffusive waiting time den- 
sities. Our analysis illustrates the inability of the FFPE to correctly predict the higher order 
moments. We demonstrate how this discrepancy creeps in while taking the continuum limit 
from the GME. For sinusoidally varying fields, the effects of the additional terms on the 
skewness and kurtosis of the position probability density were studied. We also discussed 
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how the predictions of the model can alter for different choices of the jump probability 
density and possible situations where these become relevant. 

However, whether the results correspond to physical reality or not would depend on the 
extent to which the decoupled CTRW model itself can describe natural phenomena. Possible 
areas of concern include non-ergodicity [31], dependence of the waiting time statistics on 
the external field |21], non-negligible travel time in jumps etc. Although most recent works 
concentrate on the integral moments of the position probability density, it has been suggested 
that fractional moments can give important information about its scaling behavior [32]. It 
would also be interesting to find whether the direct probabilistic approach can be extended to 
incorporate both space and time dependence of external fields. Thus the CTRW formalism, 
which has been successful in modelling many natural phenomena, still requires plenty of 
theoretical and experimental research. 

Appendix A: Calculation of A7/7 



For an oscillating external field given by, /(t) = /osinci;t, (41) and (62) give 



< a; > (t) = ^ [ (e'"" - e-'"")M^(r)rfr . (Al) 



2^ Jo 

Taking Laplace transform and applying shifting property, one obtains, 

< X >(s) = ^{Mjs - ico) - MJs + lu)} . (A2) 
2is 

For s — )■ (corresponding to t — )■ 00), this becomes, 

< X >{s) ~ -/i/o Im{M^(— iw)} or, <x> {t) ~ fifo lm{M^{-iu)} . (A3) 



For n = 2, (|64|) yields, 

M^T -r) <x> {r')dT'dT . (A4) 



< x^ > (t) = [ MJT)dT + 2fifQ [ sinwr^ / 
Jo Jo dr Jo 



Application of the shifting and convolution property of Laplace transforms on the above 
equation gives, 

< x2 >{s) = + t^Jl(^M^(^s -iu)< X >{s-iuj) - M^{s + iu)< x >{s + iu)) . (A5) 

s %s 
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Using (A3), we get, for s — )■ 0, 



< >{s] 



+ fi'f^Re{M4-tuj)} 



M4s 



- -fi'f^Re{M4-tu)M4-t2uj)} . (A6) 



Hence, for large t, 

<x^> (t) ~ {l+fi^f^Re{M4-tu)}) [ M4T)dT -fi^f^Re{M4-tu)M4-t2u)} . (A7) 

Jo 

For subdiffusive waiting time densities M^{t) ~ t"~^, which means that the first term in 



(A7) increases as and dominates the behavior of < > (t). 

Substituting n = 3 in (64), one arrives at the following equation for < 

< > it) = /i/o I s'muTM^{T)dT + 3/i/o [ sinwr— f M^{t — r') < x"^ > {T')dT'dr 



+ 3 M^it-r) <x> {T)dT 
Jo 



(A8) 



We want to analyze the effect of the first term in (A8) on the skewness. Comparison with 



(Al ) shows that this is identical to < x >. Thus, according to (A3), for large value of t, this 



approaches the constant fifo Im{M^(— za;)}. 



The 3rd term in (A8) can be written as. 



X >(s 

Jo 



3/i/o 
2i 



which, for s — > becomes 3^fQ\'m{M^{—iuj)} M^{T)dT, and increases like t° for large 
values of t. 



Let us denote the 2nd term in (AS) by a{t). In order to find out how this behaves for 



large t, we can proceed in the same way as we did in going from (A4) to (A7). This gives 



a{t) ~ 3^1^ f^lm{M^{-iuj)]Re{M^{-iuj)] / M^{T)dT + 6fifolm{M^{-iuj)}Re{M^{-iuj)} 

Jo 

- ^fi^f^lm{M4-tuj)M4-i2uj) {M4-i3uj) - M^-iu)) } . (A9) 

Taking the contribution of only the second and third term, one can write, 

< > (t) ~ 3/i/oIm{2i}(l + /iVoRe{2i}) / M^(r)rfr + 6fifolm{zi}Re{zi} 

Jo 



- folm{ziZ2iz3 - Zi)} 



(AlO) 
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where zi = M^^—iu), Z2 = M^{—i2uj) and 23 = M^^—iSu). Combining (A3), (A7) and 



(AlO), one obtains, for the 3rd central moment. 



fis = < > -3 < x'^ >< X > +2 < X >^ 

3 

~ 6fxfoRe{zi}lm{zi} - -fj.^f^lm{ziZ2izs - z^)} + f^lm{zi}Re{ziZ2} + 2^Vo (M-^i})" 

(All) 



The additional term in (A8), namely the first term, increases //a by an amount. 



A/i3 = fifQlm{zi}. Since the standard deviation remains unaffected, the relative increment 
of skewness is given by, 

A7 Afl3 



1 /^3 



Im{2;i} 



6Re{zi}Im{;2i} - |/i2/o'Im{ziZ2(^3 - ^1)} + 3^'' f^lm{z,}Re{z,Z2] + 2^Vo'(lm{^i}) 

(A12) 

Thus, the relative increment of skewness approaches a constant for large t. 



Appendix B: Calculation of /S.k./k 



For calculating the 4th central moment, we need to evaluate < x'* >, which is governed. 



according to (64), by the following equation. 



<:x'^ > {t)= [ MJT)dT + Ajj^fo [ sinwr-^ / MJt - t') < x > {T')dT'dT 
Jo Jo dr Jo 

+ 6 / M^{t -t) <x^ > {r)dT + 4^/0 / sin cur— / M^(r - r') <x^ > {T')dT'dT 
Jo Jo dr Jo 



(Bl) 



The first and second terms in (Bl) arise in our analysis but are absent in Ref. We 
have already encountered these terms in Appendix [A] while evaluating < x^ > (see (A4)). 



The results we obtained there in (A7) show that for large values of t, they increase the 4th 
central moment /i4 by an amount, 

(A/X4)i ^ (1 + 2/iVo'Re{zi}) / M^(r)dr - 2fi^f^Re{z,Z2} . (B2) 

Jo 



23 



The third term in (Bl) (call it b{t)) is just the convolution of M^{t) and < > (t) and 



hence, using (A6), one can write, 



= 6(1 + iJ.'fSRe{z,]) f M,^(t) » Mt(T)dr - f^Reiz^z.^} f M,f(T)dT . (B3) 

Jo Jo 

For subdiffusive waiting time densities, M^{T)dT ~ and M^(r) * M^{T)dT ~ t^". 

Hence h{t) ~ t^a^ 

Application of shifting and convolution property of Laplace transforms on the last term 



in (Bl) (call it c(t)) yields, after substituting results from Appendix |A 

C(t) = Co(t) + Cadd{t) , 



(B4) 



where Cadd{t) is the contribution of the additional term in < >, which increases the 4th 
central moment by (A/i4)2, given by, 

(A/i4)2 = Caddit) ^ 2fi^f^Re{zi} [ M^{T)dT - 2^? flRe{z^Z2} , (B5) 

Jo 

and, co(t) is the contribution of rest of the terms, given by, 
coit) ~ 6/iVoRe{^i}(l + /xVoRe{;2i}) / M^(r) * M^(r)dr 



+ [2Re{zi2} - /iVo'(Re{2i}Re{;2iZ2} + lm{zi}lm{ziZ2})] / M^(r)dr 

- Q^? fQRe{ziZ2{zi + 2:2)} + 3yu''/o Re {2:1 2:2 [2:3(^4 - ^2) - ^1^2]} • (B6) 

where Z4 = M^(— 240;). 

We know that the 4th central moment can be written as. 



/i4 = < X > -4 < x >< X > +6 < X >< X > -3 < X > . 



(B7) 



Consequently, the additional term in < x^ >, namely, < x > (see (A8)), directly affects 



/i4 through the second term in (B7) as, 

(A/X4)3 = -4 < X >2 ~ -4/i2/o'(Im{;2i})2 . 



(B8) 



Combining (B2), (B5) and (B8), we find the total increment of /i4 due to the extra terms 
in our analysis, 

A/i4 ^ (1 + 4/iVoRe{^i}) / M^(r)rfr - 4;uVoRe{^i^2} - 4^Vo'(M^i})' • (B9) 

Jo 
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Similarly, substituting from (A3), (A7), (AlO), (B3) and (B6) into (B7), one obtains the 
4th central moment as, 

/i4 ^ 6(1 +/iVoM^i})' / M^{t) * M4T)dT + A [ M4T)dT + B. (BIO) 

Jo Jo 

where A and B are constants related to Zi {i = 1,2,3,4). For subdiffusive waiting time 
densities, the first term increases as t^" and the second as t". Therefore the relative increment 
of kurtosis is given by, 
Ak A/i4 



/i4 



(1 + 4/i2/o Re{zi}) /o M^{T)dT + C 



(Bll) 



6(1 + /iVoM^i})' /o M^r) * M4r)dT + A /J M^(r)dr + B ' 
where C denotes —4:^'^fQRe{ziZ2} — ifi^ /q (lm{zi}Y . It is evident from this expression that 
the relative increment of kurtosis decays as t~". 
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